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<U , Abstract: 



The decay of the Standard Model Higgs boson into four leptons via a virtual W- 
boson or Z-boson pair is one of the most important decay modes in the Higgs-boson 
/S ' search at the LHC We present the complete electroweak radiative corrections of 0{a) to 

c^ • these processes, including improvements beyond 0{a) originating from heavy-Higgs effects 

and final-state radiation. The intermediate W- and Z-boson resonances are described 
(without any expansion or on-shell approximation) by consistently employing complex 
mass parameters for the gauge bosons (complex-mass scheme). The corrections to partial 
decay widths typically amount to some per cent and increase with growing Higgs mass Mh, 
reaching about 8% at Mh ~ 500 GeV. For not too large Higgs masses (Mh ^ 400 GeV) 
the corrections to the partial decay widths can be reproduced within < 2% by simple 
approximations. For angular distributions the corrections are somewhat larger and distort 
the shapes. For invariant-mass distributions of fermion pairs they can reach several tens 
of per cent depending on the treatment of photon radiation. The discussed corrections 
have been implemented in a Monte Carlo event generator called Prophecy4f.''" 
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TThe computer code can be obtained from the authors upon request. 



1 Introduction 

The primary task of the LHC will be the detection and the study of the Higgs boson. 
If it is heavier than 140 GeV and behaves as predicted by the Standard Model (SM), it 
decays predominantly into gauge-boson pairs and subsequently into four light fermions. 
From a Higgs-boson mass Mh of about 130 GeV up to the Z-boson-pair threshold 2Mz, 
the decay signature H(— > WW) — > 2 leptons + missing px [11112] has the highest discovery 
potential for the Higgs boson at the LHC [E]- For higher Higgs masses, the leading 
role is taken over by the "gold-plated" channel H — *^ ZZ — >^ 4 leptons, which will allow 
for the most accurate measurement of Mh above 130 GeV [0]. More details and recent 
developments concerning Higgs studies at the LHC can be found in the literature [El El El- 
At a future e"'"e~ linear collider [IHUHICSI; the decays H ^ 4/ will enable measurements 
of the H —>■ WW/ZZ branching ratios at the level of a few to 10% [ITT]. 

A kinematical reconstruction of the Higgs boson and of the virtual W and Z bosons 
requires the study of distributions defined from the kinematics of the decay fermions. 
Thereby, it is important to include radiative corrections, in particular real photon radia- 
tion. In addition, the verification of the spin and of the CP properties of the Higgs boson 
relies on the study of angular, energy, and invariant-mass distributions [[121 [13] • In partic- 
ular, the sensitivity of the angle between the two Z-decay planes in H ^ ZZ — ;► 4 leptons 
has been frequently emphasized in the literature. As a consequence a Monte Carlo gen- 
erator for H —>■ WW /TIL -^ 4 fermions including all relevant corrections is needed. 

The theoretical description of the decays of a SM Higgs boson into W- or Z-boson pairs 
started with lowest-order formulas for the partial decay widths. The first calculations [ 
IT^ that include off-shell effects of the gauge bosons made the approximation that one 
of the W or Z bosons was still on shell, an approximation that turns out to be not 
sufficient. Later calculations [il5j dealt with the situation of two intermediate off-shell 
gauge bosons. The various approaches are compared, e.g., in Ref. [US]- We note that 
the program Hdecay [E], which is frequently used in practice, calculates the partial 
decay widths for H — ^ WW/ZZ with on- or off-shell gauge bosons depending on Mh. 
Distributions of the decay fermions have been considered in Refs. [ ^l E] , but still in 
lowest order of perturbation theory. 

In the past the electroweak 0{a) corrections to decays into gauge bosons, H -^ 
WW /TIL, were known [^1^] only in narrow-width approximation (NWA), i.e. for on- 
shell W and Z bosons. In this case, also leading two-loop corrections enhanced by powers 
of the top-quark mass [[201 [21] or of the Higgs-boson mass [[221 [22] have been calculated. 
However, near and below the gauge-boson-pair thresholds the NWA is not applicable, so 
that only the lowest-order results exist in this Mh range. Recently electroweak corrections 
to the processes H -^ WW/ZZ -^ 4/ with off-shell gauge bosons have been considered. 
Progress on a calculation of the electromagnetic corrections to H — >■ ZZ — ;> 4 leptons has 
been reported at the RADCOR05 conference by Carloni Calame [El]- There we have also 
presented first results of our calculation of the complete 0{a) corrections to the general 
H -^ 4 leptons processes [[23- 

In this paper we describe the details of our calculation of the 0{q) corrections and 
of the included improvements beyond this order. The involved Feynman diagrams are 
closely related to the ones of the production process e+e~ — > z/z/H, whose electroweak 



0{a) corrections have been evaluated in Refs. [123 123 • Therefore, we proceed in the 
algebraic reduction of the one- loop diagrams as described in Ref. [I2Z|- On the other hand, 
the resonance structure of the decays H — > WW/ZZ -^ 4/ is practically the same as in 
e"'"e~ — > WW — > 4/, which was treated at the one-loop level in Ref. [^. Thus, we apply 
the "complex- mass scheme" [ |2Hl 12^] , where gauge-boson masses are consistently treated 
as complex quantities. This procedure fully maintains gauge invariance at the price of 
having complex gauge-boson masses everywhere, i.e. also in couplings and loop integrals. 
For a numerically stable evaluation of the latter we employ the methods described in 
Refs. [ Eni 1^ • The combination of virtual and real photon corrections is performed in 
the dipole subtraction approach [ |22l IHB] and checked by the alternative of phase-space 
slicing. The whole calculation has been implemented in a Monte Carlo generator called 
Prophecy4f. 

The paper is organized as follows. In Section |21 we fix our conventions and give ex- 
plicit results for the tree- level amplitudes of the processes H — > WW/ZZ -^ 4/. Sectional 
contains a description of our calculation of the virtual one-loop corrections; the real pho- 
ton corrections are considered in Section |3J Some details on the employed Monte Carlo 
techniques are given in Section In Section IHl we construct an "improved Born ap- 
proximation" (IBA) which approximates our state-of-the-art prediction for partial decay 
widths within < 2% for not too large Higgs masses, Mh ^ 400 GeV. Our numerical re- 
sults are discussed in Section |71 comprising partial decay widths of several representative 
H — i> WW/ZZ — >■ 41 channels as well as differential cross sections for selected channels 
in invariant masses of lepton pairs and in various angles. A comparison with results ob- 
tained from Hdecay for the partial widths is also performed there. Section |S1 contains 
our conclusions. 

2 Conventions and lowest-order results 

We consider the lowest-order processes 

H(p) ^ Mh, (Ti) + f2{k2, 0-2) + Mh, a,) + Mh, CT,), (2.1) 

where the momenta and helicities of the external particles are indicated in parentheses. 
The helicities take the values ai = ±1/2, but we often use only the sign to indicate the 
helicity. The masses of the external fermions are neglected whenever possible; they are 
only taken into account in the mass-singular logarithms originating from collinear final- 
state radiation (FSR). The matrix elements can be constructed from the generic diagram 
shown in Figure ^ 

Using the conventions of Refs. [ OH El] we denote the relevant couplings in the follow- 
ing by 



9jff = -Qf, 9zff = — -Qf, 9zff = — -Qf + 
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Figure 1: Generic lowest-order diagram for H ^ 4/ where V = W, Z. 

where Qf is the relative charge of the fermion /, and I^ f = ±1/2 the third component 
of the weak isospin of the left-handed part of the fermion field /. The CKM matrix has 
been consistently set to the unit matrix, which has no sizeable effects on our results. In 
(|2.2p we have already indicated that we use complex gauge-boson masses fiy everywhere, 

f^l^ = M^- iMvTv, \/ = W, Z, (2.3) 

where My and Ty denote the real pole-mass and width parameters. Accordingly the sine 
and cosine of the weak mixing angle are fixed by 

cl = l-sl = ^, (2.4) 

i.e. Cw and Sw are complex quantities. More details about the complex-mass scheme are 
given in Section 13.2.21 

In order to express the amplitudes in a compact way, we use the Weyl-van der Waerden 
(WvdW) spinor technique as formulated in Ref. [EE]- The spinor products (...) are 
defined by 



(pq) = e paQb = 2v/p^ 



e-i'/'p cos ^ sin ^ - e-'*^' cos ^ sin ^' 



2 2 2 2 

where pa, Qa are the associated momentum spinors for the light-like momenta 



(2.5) 



P^ = Po(l) sin 9p cos (pp, sin 9p sin (pp, cos 9p) , 

1^ = qo{l, sin 9q cos (pq, sin 9q sin (pg, cos 9q). (2.6) 

In the notation of Ref. [i36| the generic lowest-order amplitude reads 

_^yy,.„.,..a,(^^^ fc,, K, k,) = 2e'g-yX0^y^f^f,guvv 6„.,-.A.,-^AlLika, h, h, k,), (2.7) 

or more specifically for the case of Z-mediated and W-mediated decays 



3 „0-a ^CTc 



^ww,........^^^^ A:,, fc„ k,) = ^ 5.,,_5.„+5<,,,_5.„+ AW_^(fc„ k,, ke, k,). (2.8) 



The auxiliary functions are given by 

A [ka, fcfc, fcc) l^d) = 



{hkd)* {kakc) 



[{ka + k,y-fiU{kc + k,y-fil. 



21' 



^-j — yf^ai kbi kci k(i) — A — [kf), ka, kci kd) , 
■^ — \-[kai kh^ kc^ kd) = A — [ka, ki), kd, kf.), 



^4.+ (,n^a; l^b, l^c, kd) — A [kb, ka, kd, k^) , 

and obey the relations 

A_^^_^^[ka, kb, kc, kd) = yA^^^^yka, kb, kc, kd) 

^-aa,ac\ka.-, kb, kc, kd) = A^^^^[kb, ka, k^, kd), 

^aa-(Tc\'^a, kb, kc, kd) = A^^^^yka, kb, kd, kc), 

^aa<Tc\ka,j kb, kc, kd) = yA^^^^ykb, ka, kd, kc) 

^(Ta(Tc\ ai kb, kc, kd) = A^^^^ykc, kd, ka, kb)- 



(2.9) 



Mv^Mv 



MV"*Mv 



(2.10) 



The relations between the A" functions that differ in all helicities result from a P trans- 
formation. Those where only one fermion helicity is reversed are related to C symmetry. 
The last but one is due to CP symmetry, and the last one results from a symmetry under 
the exchange of the two fermion pairs. The replacements fiy ~^ A^y iii ()2.10|) ensure that 
the vector-boson masses remain unaffected by complex conjugation. 

From the generic matrix element A^g '"'''^'''^''"'^ika,kb,kc,kd) the matrix elements for 
the specific processes can be constructed as follows. To write down the explicit matrix 
elements for the different final states, we denote different fermions (/ 7^ F) by / and F, 
and their weak-isospin partners by /' and F', respectively: 



H ^ ffFF: 



H -^ ff'FF': 



Mr'^'^'^ki, k2, ks, k,) = Mf^''"''''"'\ki, k2, ks, k^) 



J^aia2a-j.ai 



{ki,k2,k^,ki) 



«^WW,o-iCT2(T30-4 



{ki,k2,k'i,ki) 



(2.11) 



(2.12) 



H - ////: 



Ml'''''''''\k,,k2,k^,k^) 



Mf^''^^^''^^\k,,k2,k^,k^) 



^0 



ZZ,(Ticr4(T30"2 



{ki,ki,kz,k2). 



(2.13) 



H - fff'f: 



Ml''''''''''\h,k2,k^,k^) 



/ZZ,cri(T2iT3cr4 , 



-M'^'^'''^''''''''\k„k„k„k2) 



(2.14) 



The relative signs between contributions of the basic subamphtudes to the full matrix 
elements account for the sign changes resulting from interchanging external fermion lines. 
The matrix elements of (j2.11|) and (j2.12|) can be extended to the case of semi-leptonic 
or hadronic final states by simply multiplying the squared matrix element by a colour 
factor 3 or 9, respectively. On the other hand, care has to be taken in the cases of ()2.13p 
and ()2.14|) for hadronic final states (semi-leptonic final states do not exist) owing to the 
non-trivial colour interferences. Summing over the colour degrees of freedom, we have 

• H — i> qqqq: 
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+ 9 



Mf^-^'^^-^'^\k,,k,,h,h) 



-QRe[M 



ZZ,o-i(T20"3'T4 





{ki, k2, k^, k,) (M^,^'^'''"''''\k,, ^4, k^, k2)y} , (2.15) 



H 



qqq'q': 



\M 



0"lCT2 0-30-4 





(ki, k2, ks, ki)\' 



M 



ZZ,cri(J2(T3Cr4 




{ki,k2,k3,ki 



+ 9 



-Mo 



WW,CTl(T4(T3fT2 



(ki,k4,k3,k2 



- 6Re {A<f '"^"^"^"^(fci, k2, ks, k^) (M 



WW,(Ti(T4Cr3(T2 




(A;i,A;4,A;3,fc2))*}- (2-16) 



Here q denotes any quark of the first two generations and q' its weak-isospin partner. 

Having constructed the matrix elements, we can write the lowest-order decay width 
To as 



dr 







2Mh 



a"l,o-2,o-3,(T4=±2 

where the phase-space integral is defined by 

4 " d^k,- ^ 



d$o = n 



^i=l 



(2.,»2.oj (-)*^'n-i:^^ 



(2.18) 



For the case H -^ ffff, which involves two pairs of identical particles in the final state, 
we implicitly include a factor 1/4 in the phase-space integral, without making this factor 
explicit in the formulas. 

3 Virtual corrections 



3.1 Survey of one-loop diagrams 

The virtual corrections receive contributions from self-energy, vertex, box, and pen- 
tagon diagrams. The structural diagrams containing the generic contributions of vertex 
functions are summarized in Figure |21 Here and in the following we omit all diagrams 
that vanish in the limit of vanishing external fermion masses from the beginning. For 
charged- current processes the generic field V stands for the W-boson field, for neutral- 
current processes we have V^ = Z, 7, where the photon is absent in couplings to the Higgs 





H 




H 




H 




Figure 2: Generic contributions of different vertex functions to H — > WW/ZZ 
where tlie blobs stand for one-particle-irreducible one-loop vertex functions. 



4/, 



boson. The generic diagrams cover all structures relevant for electroweak corrections to 
arbitrary four-fermion final states, including quarks. Note, however, that some four-quark 
final states receive corrections from diagrams with intermediate gluons on tree-like lines 
(quark- loop-induced Hgg vertex). Possible QCD corrections for quarks in the final state 
will not be considered in the following lists of diagrams. 

The pentagon diagrams are shown in Figures El and El respectively. The specific subdi- 
agrams of loop- induced 4-point functions have been shown in Ref. [I^Zl? where the process 
class e"'"e^ -^ vvYi was analyzed at one loop. They involve 4-point vertex functions of 

the type uiUiZB., z//z/;7H, /~Z+ZH, /~Z+7H, and /^ u iW^B. with I = e,fi,T denoting any 
charged lepton. The diagrams for the /"/"'"7H vertex function can be obtained from those 
for the Z~/+ZH vertex function by replacing the external Z boson by a photon and omit- 
ting the diagram where the photon couples to neutrinos. The 3-point loop insertions in 
the B.uiUi, H/^/+, HWW, HZZ, and HZ7 vertices have also been listed in Ref. [I2Z|; the 
one-loop diagrams for the H77 vertex follow from the HZZ or HZ7 case by obvious sub- 
stitutions and omissions. Most of the diagrams for the self-energies and the viViL., l~l^Z, 

z/^W^ vertex functions can be found in Ref. [IHTj. 



and /^ 







Figure 3: Pentagon diagrams for H ^ ZZ ^ ffFF, where / and F are different fermions 
with respective weak-isospin partners /' and F'. The fermion arrows in the diagrams 
involving W bosons apply to the case where both / and F have weak isospin I^ f = +1/2; 
for fermions with /^ r ■- 



■1/2 the corresponding fermion arrows have to be reversed. 







Figure 4: Pentagon diagrams for H -^ WW -^ ff'FF', where / and F are different 
fermions with respective weak-isospin partners /' and F'. The fermion arrows apply to 
the case where / and F have weak isospin I^ f = +1/2 and I^ p = —1/2, respectively; for 
fermions with opposite weak isospin the corresponding fermion arrows have to be reversed. 



All pentagon and box diagrams as well as the H77 vertex function are UV finite; also 
the B.uii'i and H/"*"/" vertex functions are UV finite since we neglect the masses of the 
light fermions everywhere apart from the mass-singular logarithms. For the other vertex 

functions, HWW, HZZ, HZ7, z/^z/^Z, l~l~^Z, l^ uiW^, and for the relevant self-energies the 
corresponding counterterm diagrams have to be included. 

3.2 Calculation of the one-loop corrections 

3.2.1 Algebraic reduction of diagrams and standard matrix elements 

The algebraic part of the two calculations has been carried out in the same way as 
in the one-loop calculation of e"'"e~ —>■ uuB. described in Ref. [I2Z1- This means that we 
separate the fermion spinor chains from the rest of the amplitude by defining standard 
matrix elements (SME). To introduce a compact notation for the SME, the tensors 

^tap-y} = '^fcikc) {7a, 7a7p7l} ^rVf^{kd) (3.1) 

are defined with obvious notations for the Dirac spinors Uf^{ka), etc., and u!± = (1 ±75)72 
denote the right- and left-handed chirality projectors. Here and in the following, each 
entry in the set within curly brackets refers to a single object, i.e. from the first line in the 
equation above we have Tf'" = UfSka)la^aVj^{h) and T^^^^ = UfSkahalnl-^^aVf^ih), 
etc. Furthermore, symbols like Tp are used as shorthand for the contraction T^p^^. We 
define the 52 SME 



^{1,2} 


"pafe,o- -pcd,T,{a,akaki,} 


i''xabcd,crT 
-^{3,4} 


-pab,a -pcd,T,{a,akakt,} 


"jabcd^cFT 
^{5,6} 


^ Yab,(7 Ycd,T,{kaM 

kc ' 


i''jabcd,crT 
-^{7,8} 


_ -pab,a -pcd,T,{ka,kt} 
fed ' 


'\iabcd,ijT 
'{9,10} 


_ -pab,a -pcd,T,{af3ka,aPk,,} 
-*- af3kc '- ' 


•'^{11,12} 


_ -pab,a -pcd,T,{af3ka,aPk,,} 


"jabcd^uT 
^13 


— ■■- a/37 







(3.2) 
The SME are evaluated within the WvdW spinor technique, similar to the lowest-order 



amplitudes described in the previous section. The tree-level and one-loop amplitudes 
A^o ^ '"'^ ^^d M-l '^ '"'^ , respectively, for the generic ft 
be expanded in terms of linear combinations of SME, 



A^o ^ '"'^ ^^d M-l '^ '"'^ , respectively, for the generic four-fermion final state fafbfcfd can 



_^aicd,ar ^ ^ p.abcd,ar ^abcd,ar ^ n = 0,l, (3.3) 

with Lorentz-invariant functions F^^ '"'^ . In this notation the lowest-order amplitudes 
(EH) read 

-'^'■0 K'^a^ '^b, t^c, l^d) — ^ "o-a,-o-b<^(Tc,-o-d 

y, i/\jabcd,aacrc 

[{ka + hy-f,l][ik, + k,)^-f,l]' ' 



yV.WW,cra<T6<7cad/T , , , N _ ^JiW r ^ J) J) 

X i iCAabcd, /o 4\ 

For the one-loop amplitudes in general all invariant functions receive contributions, in 
particular, they contain the loop integrals. The one-loop amplitudes for the various final 
states are constructed from the amplitudes for H -^ ffFF and H -^ ff'FF' as described 
in fj2.11|) to ()2.14j) for the lowest order. The one-loop correction to the partial decay 
widths, finally, reads 



/ 



dr^irt = ^/d<^o E 2Re{^r''^^'^^"''(-M^^'-"^«'-)*} . (3.5) 



0"l,0"2,0"3,(T4=±:j 
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The actual calculation of the one-loop diagrams has been carried out in the 't Hooft- 
Feynman gauge. The Feynman graphs are evaluated in two completely independent ways, 
leading to two independent computer codes. The results of the two codes are in good 
numerical agreement (i.e. within more than 10 digits for non-exceptional phase-space 
points). 

In the first calculation, the Feynman graphs are generated with FeynArts version 1.0 [ 
IHSj . With the help of Mathematica routines the amplitudes are expressed in terms of SME 
and coefficients of tensor integrals. The output is processed into a Fortran program for the 
numerical evaluation. This calculation of the virtual corrections has been repeated using 
the background- field method [^^, where the individual contributions from self-energy, 
vertex, and box corrections differ from their counterparts in the conventional formalism. 
The total one-loop corrections of the conventional and of the background-field approach 
were found to be in perfect numerical agreement. 

The second calculation has been made using FeynArts version 3 [ El] for the gener- 
ation and FormCalc [HU] for the evaluation of the amplitudes. The analytical results of 
FormCalc in terms of Weyl-spinor chains and their coefficients have been translated to 
C++ code for the numerical evaluation. 

3.2.2 Gauge-boson resonances and complex-mass scheme 

The description of resonances in (standard) perturbation theory requires a Dyson 
summation of self-energy insertions in the resonant propagator in order to introduce 
the imaginary part provided by the finite decay width into the propagator denominator. 
This procedure in general violates gauge invariance, i.e. destroys Slavnov-Taylor or Ward 
identities and disturbs the cancellation of gauge-parameter dependences, because different 
perturbative orders are mixed (see, for instance, Ref. [E] and references therein). 

In both of our two calculations we employ the so-called "complex-mass scheme" , which 
was introduced in Ref. [EH] for lowest-order calculations and generalized to the one- loop 
level in Ref. [I2H]- In this approach the W- and Z-boson masses are consistently con- 
sidered as complex quantities, defined as the locations of the propagator poles in the 
complex plane. To this end, bare real masses are split into complex renormalized masses 
and complex counterterms. Since the bare Lagrangian is not changed, double counting 



does not occur. Perturbative calculations can be performed as usual, only parameters 
and counterterms, in particular the electroweak mixing angle defined from the ratio of 
the W- and Z-boson masses, become complex. Since we only perform an analytic contin- 
uation of the parameters, all relations that follow from gauge invariance, such as Ward 
identities, remain valid. As a consequence the amplitudes are gauge independent, and 
unitarity cancellations are respected. Moreover, the on-shell renormalization scheme can 
straightforwardly be transferred to the complex- mass scheme [I28j. 

The use of complex gauge-boson masses necessitates the consistent use of these com- 
plex masses also in loop integrals. The scalar master integrals are evaluated for complex 
masses using the methods and results of Refs. [ |121 HSl HI] • 

We also treat the width of the top quark in the complex-mass scheme, i.e. we introduce 
a complex top-quark mass /xt via /i^ = in^ — imtFt. 

3.2.3 Numerically stable evaluation of one-loop integrals 

The one-loop calculation of the decay H — > 4/ requires the evaluation of 5-point one- 
loop tensor integrals. We calculate the 5-point integrals by directly reducing them to five 
4-point functions, as described in Refs. [E3 E]- Note that this reduction does not involve 
inverse Gram determinants composed of external momenta, which naturally occur in the 
Passarino-Veltman reduction [|^ of tensor to scalar integrals. The latter procedure leads 
to serious numerical problems when the Gram determinants become small. 

Tensor 4-point and 3-point integrals are reduced to scalar integrals with the Passarino- 
Veltman algorithm [|^ as long as no small Gram determinant appears in the reduction. If 
small Gram determinants occur, two alternative schemes are applied [^^. In one method, 
we evaluate a specific tensor coefficient, the integrand of which is logarithmic in Feynman 
parametrization, by numerical integration. Then the remaining coefficients as well as 
the standard scalar integral are algebraically derived from this coefficient. This method 
is used in the first loop calculation described in Section [3.2.11 The second, alternative 
method, which is used in the second loop calculation described in Section [3.2.11 makes use 
of expansions of the tensor coefficients about the limit of vanishing Gram determinants 
and possibly other kinematical determinants. In this way, all tensor coefficients can be 
expressed in terms of the standard scalar functions. 

The whole procedure for the evaluation of the scalar and tensor one-loop integrals has 
been taken over from the one- loop calculation of e"'"e~ -^ 4fermions [I28j. 

3.2.4 Input-parameter scheme 

We use the "G^ scheme", where a large universal part of the 0{a) corrections is 
absorbed into the lowest-order prediction. In this scheme the electromagnetic coupling 
constant a = e^/(47r) is derived from the Fermi constant G^, the muon decay constant, 
according to 



n Ml 



«c. = ^^^^^£^ ( ^ - "^ j • (3-^) 



This procedure takes into account the running of the electromagnetic coupling constant 
a(Q^) from Q^ = to the electroweak scale and also accounts for universal corrections 
related to the p parameter in the coupling of the W boson to fermions. 
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In order to avoid double-counting, the corrections absorbed in the lowest-order predic- 
tion by the use of ac^ have to be subtracted from the explicit 0{a) corrections. To this 
end, we subtract the weak corrections to niuon decay Ar [ EH US] from the corrections 
in the a(0) (on-shell) scheme. This can be done by redefining the charge renormalization 
constant as 



6Z, 



'g„ 



6Z, 



.(o)-^(^^)— ' (3.7) 



where (Ar)i_ioop is the one-loop expression for Ar evaluated in the complex- mass scheme. 

3.3 Leading two-loop corrections 

Since corrections due to the self-interaction of the Higgs boson become important for 
large Higgs masses, we have included the dominant two-loop corrections to the decay 
H —>■ VV in the large-Higgs-mass limit which were calculated in Refs. [122 123 • They are 
of order 0{GlM^) and read 

j dr^.M^ = 62.0308(86) ( j^ 1 / dlo, (3.8) 



where the numerical prefactor has been taken from Ref . [ |22] • The error of this factor is 
far beyond other uncertainties and, thus, ignored in the numerics. 

We do not include any higher-order corrections proportional to a power of G^m^ since 
we already see at the one-loop level that the heavy-top limit does not provide a sound 
approximation of the corrections from closed fermion loops. In particular, for Higgs masses 
near and above the tt threshold, Mh ^ 2mt the large-mt limit is not appropriate. 

4 Real photon corrections 

4.1 Matrix element for H — > 4fj 

The real photon corrections are induced by the process 

H(p) — > fiikuai) + f2{k2,a2) + f3{h,(r3) + h{h,cT4) +-f{k,X), (4.1) 

where the momenta and helicities of the external particles are indicated in parentheses. 

As for the lowest-order process, we consistently neglect fermion masses. However, we 
restore the mass-singular logarithms appearing in coUinear photon emission as described 
in Section W?^ 

The matrix elements for the radiative process can be constructed in the same way as 
for the lowest-order process ()2.H) from the set of generic diagrams that is obtained from 
Figure Q] by adding a photon line in all possible ways to the charged particles (including 
possible new graphs involving would-be Goldstone-boson exchange). 

We have evaluated the generic helicity matrix elements A4'^'''^'''^'''^'^^{ka, kb, kc, kd, k) of 
this process again using the WvdW spinor technique in the formulation of Ref. [ Q^ . The 
amplitudes generically read 

^VV,a.a,a.a,X^Q^^ Qfc, Qc, Qd, ka, h, k,, k^, k) = (4.2) 

2^26^ 9Vfah9v%h9^VV Sa^-at^a.-aa ^"^^a.xiQa, Qb, Qc, Qd, ka, h, k^, kd, k) , 
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or more specifically for the case of Z-mediated and W-mediated decays 

A1^ ' " " [Qa^ Qbj Qci Qdi ^aj "^65 ^c; ^dj "^j = 






5„ 



-o-i,0(j^^_(jj 



^o-aCTcAl^a? Qbi Qci Qdi ka, Kb, Kc, Kd, k), 



■^y ' " ^ [Qa, Qb) Vc) Qdi ka, fcfe, fcc? "^d) k) 



V^e'^/iw 



Oo-a,-"o-i,,+0o-c, — "o-d,+ ^ X\Qa, Qb, Qc, Qd, ka, kb, Kc, kd, k). 



(4.3) 



The auxiliary functions are given by 

^ \,Qa} Qbj Qc} Qdi ka, kb, kc, kd, k) = 

{kahYikakc) + {kkbY{kkc) 



{kbkdY 



[{ka + kb + kf- fll] [{kc + kdf- fil] 



X 



Qa 



+ 



Qa — Qb 



{kka) 



{kka)*{kkb)* {ka + kb)^ - fil {kkb)* 

{kckdY{kcka) + {kkdY{kka) 

[{ka + kbY-^ll][{kc + kd+kf-^ll] 

Qc , Qc- Qd {kkc) 



X 



{kkcYikkdY {kc + kdY - fil {kkdY J 
Qa-Qb {kbkdY {kakc) 



\{ka + kbf - 111] \{kc + kdf - 14] {kkbYikkdY 



^^ [Qai Qbi Qci Qdi ka, kb, Kc, kd, k) = A [~Qbi ^Qai Qci Qdi kbi ka, kc, kd, k), 

A — I — \Qa, Qb, Qc, Qd, ka, kb, kc, kd, k) = A [Qa, Qb, Qdi Qci kai kb, kd, kc, k), 

^+^ — \Qai Qbi Qci Qdi kai kbi kci kdi k) = A {~Qbi ^Qai ~Qdi ~Qci kbi kai kdi kci k), 

o-aCTc+ W«' Qbi Qci Qdi kai kbi kci kdi k) = 

\^-aa-crc-\Qai Qbi Qci Qdi kai kbi kci kdi k) 



fJ-V^f^Y 



(4.4) 



and obey the relations 



A 



vv 

\VV 



[Qai Qbi Qci Qdi kai kbi kci kdi k) — 
^(Ta(Tc>-\^ai Qbi Qci Qdi kai kbi kci kdi k) 



Mv^Mv 



^-aa,(Tc,X\Qai Qbi Qci Qdi kai kbi kci kdi k) — 

^aa(TcX\~Qbi "Qai Qci Qdi kbi kai kci kdi k)i 

^(7a,-(Tc,X\Qai Qbi Qci Qdi kai kbi kci kdi k) = 

^aaCTcXiQai Qbi ~Qdi ~Qci kai kbi kdi kci k)i 

^o-ajO-c-ACya; Qbi Qci Qdi kai kbi kci kdi k) = 

~\^crccraX[Qdi Qci Qbi Qai kdi kci kbi kai k) 

^o-aO-cACya) Qbi Qci Qdi kai kbi kci kdi k) = 

^UraaXiQci Qdi Qai Qbi kci kdi kai kbi k) ■ 



IJ-V-^tJ'Y 



(4.5) 
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The relations between the A" functions that differ in all helicities result from a P trans- 
formation. Those, where only one fermion helicity is reversed are related to C symmetry. 
The last but one is due to CP symmetry, and the last one results from a symmetry under 
the exchange of the two fermion pairs. The charges of the fermions are related by 

Qa-Qb + Qc-Qd = 0. (4.6) 

For the Z-mediated decays, where Qa = Qb and Qc = Qd, the auxiliary function ()4.4p 
simplifies to 

-^ (.Vci) Qai Vc) Vc; ^aj ^fej ^c; "^dj "^j ^ 

{KhY{kaK) + {kkb)*{kk,) Qa 



{hkdY 



[{ka + k, + k)^ - fiUikc + ka? - fil.] {kkaY{kk,Y 

{k,kaY{kcka) + {kkaY{kka) Q, 

[{ka + hf - iil][{k, + kd + k)^ - /if.] {kk,Y{kkdY 



(4.7) 



From the generic matrix element M.}^^''^°-'^'''^'''^''^{ka, kb, kc, kd, k) the matrix elements 
for the specific processes can be constructed in complete analogy to the process without 
photon as in (J2IIH)-(I2II1I)- 

The squares of the matrix elements ()4.2j) have been successfully checked against the 
result obtained with the package Madgraph [H?] numerically. 

The contribution F^ of the radiative decay to the total decay width is given by 

^^WH J ,1 ^_ , , 

a"l,0-2,0-3,0"4 = ±2 '^— ^ti 

where the phase-space integral is defined by 

h^^-Ij^ (n / jSkf) ''^'''" i'-'-t '') ■ <"' 

Without introducing soft and collinear regulators this phase-space integral diverges in the 
soft (fc° —>■ 0) and collinear {kki -^ 0) regions. The applied solutions are outlined in the 
next section. 

4.2 Treatment of soft and collinear divergences 

In the combination of virtual and real photon corrections, the fermion-mass effects 
have to be restored in the phase-space regions of collinear photon radiation off charged 
fermions, and the IR regularization for soft-photon emission has to be implemented. To 
this end, we employ the dipole subtraction formalism for photon radiation [E2] as well as 
the more conventional phase-space slicing approach. 

The actual calculation exactly follows the one described in Ref. [ESI, where elec- 
troweak corrections to the related process 77 -^ WW -^ 4/ have been calculated. The 
structure of soft and collinear singularities of this process is exactly the same as in the 
decay H ^ 4/ considered in this work, because both processes involve the same pattern 
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of charged particles in the initial and final states. Consequently, apart from obvious sub- 
stitutions for the fiux factors all formulas given in Section 4 of Ref . [ ES] for cross sections 
literally carry over to our decay widths. 

In the matching of real and virtual corrections the issue of coUinear safety of observ- 
ables is crucial. We speak of collinear-safe observables if a nearly coUinear system of a 
charged fermion and a photon is treated inclusively, i.e. if phase-space selection cuts (or 
histogram bins of distributions) depend only on the sum ki + k oi the nearly coUinear 
fermion and photon momenta. In this case the energy fraction 

(4.10) 



kf + fco 



of a charged fermion /j after emitting a photon in a sufficiently small cone around its 
direction of fiight is fully integrated over, because it is not constrained by any phase-space 
cut (or histogram bin selection in distributions). Thus, the Kinoshita-Lee-Nauenberg 
(KLN) theorem [HH] guarantees that all singularities connected with FSR cancel between 
the virtual and real corrections, even though they are defined on different phase spaces. 
The full phase-space integration in H — >■ 4 f {+'-/), which leads to partial decay widths, 
is trivially coUinear safe. More generally, a sufficient inclusiveness is achieved by the 
photon recombination described in Section Fm which treats outgoing charged fermions 
and photons as one quasi-particle if they are very close in angle. The original version of the 
dipole subtraction formalism for photon radiation [|32] deals with collinear-safe situations 
only. The generalization to the non-coUinear-safe FSR is described in Ref. [ [33] . 

Figures and El illustrate the agreement between the subtraction and slicing methods 
for the partial decay widths of the two decay channels H -^ VeO^ ^^v^ and H — >■ e~e"'"/i~/i"'". 
These results were obtained with 5-10^ events. In the slicing approach, the phase- 
space regions of soft or coUinear photons are defined by the auxiliary cutoff parameters 
^E <^ Fw and A^ -C 1. The region k^ < AE is treated in soft-photon approximation, the 
regions 6'^/ < A^, fc° > AE {6^f is the emission angle from any fermion /) are evaluated 
using coUinear factorization. In the remaining regular phase space no regulators (photon 
and fermion masses) are used. Therefore, the slicing result is correct up to terms of 
0{AE) and 0{A6). For decreasing auxiliary parameters AE and A^, the slicing result 
reaches a plateau, as it should be, until the increasing statistical errors become large and 
eventually underestimated. In the plateau region the slicing and subtraction results are 
compatible within statistical errors, but the subtraction result shows smaller integration 
errors although the same number of events is used. 

4.3 Higher-order final-state radiation 

Photons that are emitted coUinear off a charged fermion give rise to corrections that 
are enhanced by large logarithms of the form a log m'j/M'^, where rrif is a fermion mass 
and M is some typical energy scale. If the photons are treated fully inclusively, as it is 
the case if the photons are recombined with the corresponding fermion, these logarithms 
cancel due to the KLN theorem [^1^. If, however, distributions like in the invariant 
mass of two fermions, as discussed in Section [7J are to be considered without recombining 
coUinear photons, then these logarithms do not cancel and yield large effects. Thus, 
corrections of this origin should be taken into account beyond 0{a). This can be achieved 
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Figure 5: Dependence of the relative corrections 5 to the partial decay width on the 
energy cutoff /S.E (l.h.s.) and on the angular cutoff A9 (r.h.s.) in the slicing approach for 
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obtained with the dipole subtraction method is shown as a la band in the plots. 
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in the structure-function approach [|^ which is based on the mass-factorization theorem. 
According to this theorem the decay width including the leading-logarithmic FSR terms 
can be written as 



/ 



dr 



LLFSR 



n 



dz,;rt.Lf2. M^) 



uo 



dro6cut({2;j/^i})- 



(4.11) 



The function ©cutd-Zj^j}) generically denotes all histogram routines or phase-space cuts. 
It depends on the fermion momenta Zjkj which, in the case of charged fermions, may 
be reduced by the factor Zj due to collinear photon emission. For neutral fermions we 
have Zj = 1. The structure functions including terms up to 0{a^), improved by the 
exponentiation of the soft-photon parts, read [ t5Dj 



rJ:.^''="pr2,M2) 



exp(-|A7E + |A) A^ 

r (i + IP- 

_Pl(l + 3z^ 
32 

(3f 



Z] 2 



§(1 + .) 



\\i{z) + 4(1 + z) ln(l - z) + 'b + z 



384 
+ 



[l + z) 6Li2(z) + 121n^(l -^) -37r^ 



1 + 8z + 3^2) in(^) + 6(z + 5)(1 - z) ln(l - z) 



+ 12(1 + z^) \n{z) ln(l - z) ~ -(1 + Iz^) \n^{z) 

+ -(39-24^-15^2) I (^4J2) 

with 7e and T{y) denoting Euler's constant and the Gamma function, respectively. The 
mass-singular logarithm 



A 



2Q?a(0) 



TT 



In 



\mf J 



(4.13) 



involves a scale M, which is not fixed in leading logarithmic order and should be set to 
a scale typical for the process under consideration. We use M^ = M^ in our evaluations. 
The appropriate coupling in the leading-logarithmic terms is a(0), since these originate 
from real or virtual soft or collinear photons. As the function (1 — z)^~^ is difficult to 
integrate numerically, an appropriate mapping has to be chosen. 

In order to study the influence of the higher-order terms we alternatively expand the 
exponential up to terms of 0(q;^), yielding 



r^,"(^,M2) = 5(i-^) + 



A 1 + ^' 



4 l-z 



Pf (l + Az + z' l + 3z\ . . l + z\ .^ 
32 I 1 — z V — z \ — z 



+ 



01 rl5 + 242 + 152^ 
1536 1 \-z 



47r' 



,1 + 3^2 



1 + 82; + 3^2 
o mi^z) 
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+ 24 ^ ln(l -z)- 24(1 + z) Li^iz) + 48=--— ln2(l - z) 



1-z ' ' ' / -V / ^__2 

(4.14) 



1 -I- 7^ 1 _|_ 7^2 

- 48^^— ln(z) ln(l - 2) + 2 ln'(z 

1 — z 1 — 2; 

where the [...]+ prescription is defined as usual, 

\z\f{z)\ g{z)^ f'dzf{z)[g{z)-gil)], (4.15) 

L J+ Jo 

We convolute the lowest-order width according to ()4.1ip and add this to the result for 
the (9(a)-corrected width. In order to avoid double counting, we have to subtract 



/ drLLFSR,! = / dlo + / dFo ^ / dzi Tll^'^{zi, M^)Qc^t{ziki, {kj-^i}) 



(4.16) 



i.e. the leading logarithmic terms up to 0{a), from /dFLLFSR- They are defined by 

r-'\.,M^) = ^(i±^) . (4.17) 

Note that we have to subtract the 0{a) terms according to the scheme that is applied for 
the virtual corrections. Since we work in the G^ scheme, f3i^G is proportional to ac , as 
defined in Section [7| 

5 Final prediction and Monte Carlo integration 

Summarizing all contributions to the differential decay width, ()2.17p . ()3.5|) . ()3.8p . ()4.8|) . 
fl4.11|) . and ()4.16p . we get the following prediction, 

JdT = JdTo + J dr^irt + J dTalM^ + J dT^ + J dFLLpSR - / dlLLFSR,!- (5-1) 

The phase-space integrations are performed using the multi-channel Monte Carlo tech- 
nique [ 1^ where the integrand is flattened by choosing appropriate mappings of the 
pseudo-random numbers into the momenta of the outgoing particles. In more detail, the 
Monte Carlo part of Prophecy4f builds upon the existing generators RacoonWW [ 
123 E21 and C0FFER77 [|2niES]- The results obtained with the multi-channel technique 
have been checked against a second integration program based on the adaptive integra- 
tion program VEGAS [ |H1] . Although Vegas already performs some adaption to the 
integrand structure a further improvement of the numerical error has been achieved by 
using mappings of the integration variables designed to flatten the resonances. 

The numerical results presented below have been obtained using 5-10^ events except for 
the plots showing the decay width as a function of the Higgs mass which were calculated 
using 2 ■ 10^ events per point. Since the virtual corrections (rendered flnite by adding the 
soft and coUinear singularities from the real corrections), and also their statistical error, 
are at least a factor 10 smaller than the lowest-order values for moderate Higgs masses, 
we only evaluated the virtual corrections every 100th time, which improves the run-time 
of the program but does not deteriorate the overall statistical error. Soft and collinear 
singularities were treated with the subtraction method in the results shown below. 
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6 Improved-Born Approximation 

The electroweak corrections contain large contributions of universal origin. Besides 
final-state radiation, which is discussed in Section W^ these consist, in particular, of the 
corrections associated with the running of a, corrections proportional to m^/M^, and 
corrections proportional to M^/M^. By suitable parametrization of the lowest-order 
matrix elements, some of these universal corrections can be incorporated in the lowest 
order, thus reducing the remaining corrections. This does not only reduce the 0{a) 
corrections but in general also the higher-order corrections. Corrections associated with 
the running of a, and corrections related to the p parameter in the W-boson-fermion 
coupling are incorporated in the lowest-order prediction by using the G^ scheme. 

Some loop diagrams involving top quarks lead to corrections that are enhanced by 
a large coupling factor G^m^ in the limit of a large top-quark mass rrit. For the HZZ 
and HWW vertices this type of corrections was considered even up to two-loop order in 
Ref. [E|. However, the pure heavy-top limit mt — > oo is only applicable for Mh -C 2mt. 
Since we are also interested in Mh values near and above the tt threshold, we instead 
consider the more general limit mt,MH ^ Mw,Mz,m/^t, i-e. we do not assume any 
hierarchy between the Higgs and the top-quark masses. We evaluate all closed fermion 
loops in this limit and keep only contributions that are enhanced by m^ times any func- 
tion of the ratio Mn/mt. The numerical analysis shows that this procedure yields a 
very good approximation for the sum of all closed fermion loops in the G^ scheme. In 
order to approximate the remaining bosonic corrections, the leading one- and two-loop 
corrections to the Hl^l^ vertices in the large-Higgs-mass limit, which are proportional 
to Gfj^M^ and G'tM^, respectively, are included in the IBA; these corrections are taken 
from Refs. [ |22 I2B] • Moreover, for H -^ WW -^ 4/ we include the leading effect (5coui 
of the Coulomb singularity as calculated in Ref. [ E3] , which originates from soft-photon 
exchange between the two slowly moving W bosons near the WW threshold. Finally, we 
take into account the QCD corrections to the gauge-boson decays if quarks are involved 
in the final states. The remaining corrections are expected to be widely independent of 
Higgs mass Mh and of the choice of the 4/ final state. 

Our IBA for the partial decay widths, which is constructed according to these lines, 
reads 

^ 0"l,o-2,0'3,(T4=± 



Sv^ir^ 3„ \9?,,,, g?,,,J ""Vrf 



2 



+ ^ f 5£! _ 3V3, + 15U 62.0308(86) ( ''^"'» 



2 



8v^7r2 16 2 ■ y J IiGtt^v^, 



I rQCD I rQCD I ^ 

+ ^Z^fj2 + ^Z^/3/4 + CHZZ 



/ 






-H-WW-4/ _ _J_ I ^^ I A^WW,-+- 
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2 



X ReU + 



87r2v^ 



'M 



2 ^ 

—5 + thww I — 2" 



+ V sVStt + — + 62.0308 86 _^^ 



V2. 
+ ^(^)<^coui (Ml, {h + hf, ih + hY) 

+ '^W^/i/2 + ^W^faU "^ "^HWW p (6.1) 

where the terms proportional to a charge factor Qf are absent if / is a neutrino. The 
phase-space integral was defined in ()2.18|) . The auxiliary functions rnyy, which appear 
in ()6.ip . are given by 

THZZ [ ^ j = 20 + 6A' + 3A(A' + 1) In(xt) + 3(1 - A") In'(xt), 



^h\ _« , ,o^2 , o^.o^2 ,^,.,.^ , 3 



THWW (^1=8 + 12A' + 3A(3A' - 1) In(xt) + ^{^ - A ) Inl^^t), (6.2 



where 



A 



4/it' A-1 .^ ^. 



V ^h' A + 1' 

They have the property rHvy(O) = 0, i.e. they quantify the deviation of the m^/M^- 
enhanced corrections from the pure heavy-top limit, which are made explicit in ()6.ip . 
Note that we consistently use the complex top mass /it instead of rrit. The correction 
factor ^coui containing the Coulomb singularity reads 



i^couKs, fc+, k_) = —^ Im <; In ' 



(3 I \p + p + /\ 



M , 



Js^ + ki + kt- 2skl - 2sA;2 - 2A;2_A:2 
s 



f3=Jl- ^, Am = ^^^, (6.4) 

V s s 

with the fine-structure constant a{0). The auxiliary function 

^(A = (l-^')' (6.5) 

restricts the impact of ^coui to the WW threshold region where it is valid. The QCD 
corrections 6^ , r to the gauge-boson decays are given by 



a, 



^,. = 0. ^?2„-, = ^. (6.6) 



according to whether fifj is a lepton {Ijlj) or a quark (gj^j) pair, respectively. Finally, 
the "constants" chw have been introduced to account for a sizeable, but widely Mh- 
independent offset in the relative corrections to the YiVV vertices that is induced by 
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non-leading corrections. In practice, it is often sufficient to set cuvv to a numerical 
constant, which can be determined from a comparison with the full correction to the 
H — > VV -^ 4/ process. For the input parameters given in Section mi we find that the 
choice 

chzz = 3%, chww = 4% (6.7) 

is appropriate. The values of cuvv will certainly not change significantly if the input 
parameters of Section [73 vary within their experimental uncertainties. 

When defining the IBA for the final states ffff and fff'f we make use of the fact 
that the two subamplitudes in (j2.13p and (j2.14p have a very small interference. Therefore, 
we define the IBA for the squared matrix element from the IBA of the corresponding 
squared subamplitudes and take into account the interference in lowest order without 
modification. 

In deriving the contributions enhanced by m^ in (6.1), we used a double-pole ap- 
proximation and implicitly assumed an integration over all decay angles. Without this 
assumption, corrections to angular correlations are induced by the interference of the two 
relevant formfactors describing the HV^V vertices. Hence, the IBA (6.1) is in ffist place 
constructed for partial widths. 

7 Numerical results 

7.1 Input parameters and setup 

We use the following set of input parameters [ EHl , 

6"^= 1.16637 X 10-5 GeV~^ a(0) = 1/137.03599911, as =0.1187, 

M^^ = 80.425 GeV, T^^ = 2.124 GeV, 

M^^P = 91.1876 GeV, T^^^ = 2.4952 GeV, 

me = 0.51099892 MeV, m^ = 105.658369 MeV, m^ = 1.77699 GeV, 

mu=66MeV, mc=1.2GeV, mt= 174.3 GeV, 

md=66MeV, ms=150MeV, mb = 4.3 GeV. (7.1) 

The masses of the light quarks are adjusted to reproduce the hadronic contribution to 
the photonic vacuum polarization of Ref. [E^. As discussed in Section [3.2.41 we use the 
Gn scheme, i.e. we derive the electromagnetic coupling constant from the Fermi constant 
according to (j3.6|) . We use a^,, everywhere except for the couplings of the collinear 
photons, as described in Section 1^^ In this case we use a(0), because this reflects the 
coupling behaviour of real photons. 

Using the complex-mass scheme, we employ a fixed width in the resonant W- and Z- 
boson propagators in contrast to the approach used at LEP to fit the W and Z resonances, 
where running widths are taken. Therefore, we have to convert the "on-shell" values of 
M^^P and F^^^ {V = W,Z), resulting from LEP, to the "pole values" denoted by My 
and Fy in this paper. The relation between the two sets of values is given by [ EH] 



My = M^^V^l + (FLEP/^LEP)2^ Ty = T^,''^ / s/i + Jt^^Jm^^ (7.2) 
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leading to 

Mw=80.397GeV, Mz = 91.1535 GeV. (7.3) 

We make use of these mass parameters in the numerics discussed below, although the 
difference between using My or My^^ would be hardly visible. The widths of the gauge 
bosons W and Z, Fw and Fz, are calculated from the above input including C(a) correc- 
tions, but using real mass parameters everywhere. Alternatively, the experimental widths 
calculated from (j7.2|) could be used, but the procedure pursued here ensures that the "ef- 
fective branching ratios" of the W's and Z's, which result from the integration over their 
decays, add up to one if all decay channels are summed over. The gauge-boson widths 
depend on the Higgs mass only weakly. For the Higgs masses Mh = 140, 170, 200 GeV the 
corresponding values are given in Table H These values are used everywhere in the com- 
plex masses, i.e. we also apply the 0(a)-corrected W and Z widths for the lowest-order 
predictions. 

In order to improve the corrections for partial Higgs decay widths for a Higgs mass 
near the tt threshold (Mh ~ 2mt), we evaluate loop diagrams with internal top quarks 



with a complex top-quark mass /it = \lrnl — imtFt. To this end, we set the top-quark 
width to its lowest-order prediction in the SM, 

8nV2mf 

where the simplifications mb = and Vtb = 1 are used. 

The angular distributions in Section 17.51 are defined in the rest frame of the Higgs 
boson. All observables are calculated without applying phase-space cuts, and, if not 
stated otherwise, a photon recombination is performed. More precisely, if the invariant 
mass of a photon and a charged fermion is smaller than 5 GeV, the photon momentum is 
added to the fermion momentum in the histograms. If this condition applies to more than 
one fermion the photon is recombined with the fermion that yields the smallest invariant 
mass for the resulting fermion-photon pair. 

All but the lowest-order predictions contain the higher-order FSR as described in 
Section l^ni as well as the two-loop corrections proportional to G'tM^ given in Section IT^ 

7.2 Results for partial decay widths 

In Tabled the partial decay width including 0{a) corrections is shown for different 
decay channels and different values of the Higgs mass. In parentheses the statistical 
error of the phase-space integration is indicated, and S = F/Fq — 1 labels the relative 
corrections. The first two channels, e~e~^fi~fi'^ and l^l^l^l^, I = e,iJ,, result from the 
decay H — i> ZZ — >■ 4/. The partial widths for only electrons or muons in the final 
state are equal in the limit of vanishing external fermion masses, since for collinear- 
safe observables, such as the partial widths, the fermion-mass logarithms cancel. The 
corresponding lowest-order matrix elements are given in ()2.11|) and ()2.13|) . respectively. 
The width for H -^ l~l^l~l~^ is typically smaller by a factor 2, because it gets a factor 1/4 
for identical particles in the final state and it proceeds in lowest order via two Feynman 
diagrams that are related by the exchange of two outgoing electrons and that have only a 
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Table 1: Partial decay widths for H -^ 41eptons including 0{q) and (9(G^M^) correc- 
tions and corresponding relative corrections for various decay channels and different Higgs 
masses. 

small interference. The channel z/ee"^/i~z/^ fl2.12|) results from the decay H -^ WW — > 4/, 
while the last channel vil+l-vi ()2.14|) receives contributions from both the decay into W 
and into Z bosons. The larger the Higgs mass, the larger is the decay width, because the 
available phase space grows. 

In the two upper plots of Figure [7| we show the partial decay width for the final state 
u^e+fi'u^ as a function of the Higgs mass. The lower plots show the corrections relative 
to the lowest-order result. As already explained, we always normalize to the lowest-order 
result that already includes the 0(a;)-corrected gauge-boson width in the complex masses 
of the gauge bosons. A large fraction of the 0{a) corrections is transfered to the lowest- 
order decay width by applying the G^ scheme. Thus, the corrections are at the order 
of 2-8% for moderate Higgs masses. However, for large Higgs masses the corrections 
become larger and reach about 13% at Mh = 700 GeV. In this region the leading two- 
loop corrections already amount to about 4%. Around the WW threshold at 160 GeV the 
Coulomb singularity, which originates from soft-photon exchange between the two slowly 
moving W bosons, is reflected in the shape of the curve. The influence of diagrams with a 
Higgs boson splitting into a virtual Z-boson pair (ZZ threshold) is visible at Mh ~ 2Mz. 
At about 2mt the tt threshold is visible. 

For stable W or Z bosons, i.e. in the limit Fy — > (l^ = W, Z), it is possible to 
define a narrow-width approximation (NWA) where the matrix elements factorize into 
the decay H -^ VV and the subsequent decay of the gauge bosons into fermions. By 
definition the NWA is only applicable above the WW or ZZ threshold. However, its 
analytical structure and evaluation is considerably simpler than in the case of the full 
decay H —>■ WW/ZZ -^ 4/ with off-shell gauge bosons. Therefore, above threshold the 
NWA allows for an economic way of calculating relative 0{a) corrections to the integrated 
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Figure 7: Partial decay width for H — >■ u^e'^^'u^ as a function of the Higgs mass. The 
upper plots show the absolute prediction including 0{a) and O(G^M^) corrections, and 
the lower plots show the comparison of the corresponding relative corrections with the 
NWA and IBA. 
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decay width, while the lowest-order contribution may, of course, still take into account 
unstable gauge bosons. Following this line of thought, we define 



pNWA 
'- 1 

-pNWA' 
J- n 



r^"^ = ro^;kA, (7.5) 





with 



and 



rNWA p ^V/l/2,l^V/3/4,l /7 «N 

J-1 =iHyy,i :^ — :^ , U-oJ 

J- V,!'- V,l 
pNWA _ -p ^Vfif2,0^Vf3f4fi fj 7N 

J-0 —^nvvfi 7; — 7; • ['■') 

The indices "0" and "1" label lowest-order and (9(a)-corrected results, respectively. The 
Higgs-mass-enhanced two-loop terms, described in Section EISl have also been included 
in T}ivv,i- In order to be consistent we again use the (9(a)-corrected total width for the 
gauge bosons in the denominators of the branching ratios in F^^^. We note that we have 
rederived all necessary 0{a) corrections entering the NWA; the hard-photon corrections 
to the decay H — >■ WW have been successfully checked against the expression given in 
Ref. [Iini- The NWA is evaluated with real gauge-boson and top-quark masses. 

A few GeV above the corresponding gauge-boson-pair threshold the NWA agrees with 
the complete corrections within 1%. Near Mr = 2Mz ~ 180 GeV the loop-induced 
ZZ threshold can be seen in the relative corrections to H — > WW -^ z/ee"'"/x~t'^ shown 
in Figure [7| In the NWA this threshold leads to a singularity visible as a sharp peak; 
in the off-shell calculation in the complex-mass scheme this singular structure is smeared 
out, because the finite Z-boson width is taken into account. Since the ZZ threshold 
corresponds to the situation where two Z bosons become on shell in the loop, the latter 
description with the singularity regularized by a finite Fz is closer to physical reality. An 
analogous situation can be seen near Mh = 2mt ~ 350 GeV for the tt threshold with 
top quarks in the loops. Again the inclusion of the top decay width Ft, as done in the 
off-shell calculation, yields the better description. In Figure [3 we show also the relative 
corrections obtained in the IBA of (jfi.lj) . The Coulomb singularity and the tt threshold 
are well described, and the full corrections are reproduced within < 2% for Higgs masses 
below 400 GeV. 

The plots in Figure |H1 show the decay width and the relative correction for the final 
state e~e"'"/i~/i~'". The corrections are between 2% and 4% for moderate Higgs masses and 
rise to more than 10% for large Higgs masses. At a Higgs mass of about 160 GeV the 
influence of the WW threshold can be observed. As explained above, the behaviour of 
the corrections as a function of the Higgs mass is smooth, because the finite W-boson 
width is also used in the loop integrals. In contrast to the decay H — > z/ee"'"/i~z/^, there is 
no Coulomb singularity in this channel because the Z boson is electrically neutral. The 
NWA reproduces the complete result within 0.5% not too close to the threshold. The IBA 
agrees with the complete calculation to better than 2% for not too large Higgs masses. 

7.3 Comparison of partial widths with Hdecay 

Predictions for the partial decay widths of the Higgs boson can also be obtained with 
various program packages, such as Hdecay [EI- Hdecay contains the lowest-order 
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Figure 8: Partial decay width for H -^ e~e"'"/i~/x"'" as a function of the Higgs mass. The 
upper plots show the absolute prediction including 0{a) and O(G^M^) corrections, and 
the lower plots show the comparison of the corresponding relative corrections with the 
NWA and IBA. 
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decay width for H — ^ V'^*^V^*^ and the leading one- loop corrections oc G^M^ and two- 
loop corrections oc Gf^M^. In order to obtain the decay width for H -^ WW/ZZ -^ 4/, 
we define 

J- =^HVV—F^ r ' ^'-^^ 

-L V,l J- V,l 

where r™y is the decay width from Hdecay. In ()7.8p the branching ratios of the gauge 
bosons are normalized in the same way (lowest order in the numerator, corrected total 
width in the denominator) as the effective branching ratios of our lowest-order predictions 
for the H -^ VV -^ 4/ partial widths; otherwise a comparison would not be conclusive. 

The comparison in Figure IHl where F^^ is shown relative to our complete lowest-order 
prediction, shows that Hdecay agrees with our lowest-order prediction below the decay 
threshold quite well. In this region F™y consistently takes into account the off-shell 
effects of the gauge bosons. Above the threshold Hdecay neglects off-shell effects of 
the gauge bosons. For large Mh it follows our corrected result within a few per cent, 
because the dominant radiative corrections oc G^M^ and oc G'^M^, which grow fast 
with increasing Mh, are included in both calculations. In the threshold region, off-shell 
effects are, however, very important. Here, the difference between the complete off-shell 
result and the Higgs width for on-shell gauge bosons amounts to more than 10%. In detail, 
Hdecay interpolates between the off-shell and on-shell results within a window of ±2 GeV 
around threshold. The maxima in the Hdecay curves near the WW and ZZ thresholds in 
the upper and lower left plots of Figure 121 respectively, are artifacts originating from the 
on-shell phase space of the W or Z bosons above threshold. Approaching the threshold 
from above, the on-shell phase space, and thus the corresponding partial decay width, 
tends to zero. This feature is avoided in Hdecay by the interpolation. The described 
maxima in the Hdecay curves have nothing to do with the maximum of the correction 
near the WW threshold in the upper left plot, which is due to the Coulomb singularity. 

7.4 Invariant-mass distributions 

In Figure El we study the invariant-mass distribution of the fermion pairs resulting 
from the decay of the W bosons in the decay H — > z/ee"^/i~z/^. The plots on the l.h.s. show 
the distribution for fi~i^^ including corrections for Mh = 170 GeV and Mh = 140 GeV, 
i.e. for one value of Mh above and one below the WW threshold. The plots on the 
r.h.s. compare the relative corrections for u^e'^ and /i"!^^ both with and without photon 
recombination. The invariant mass Mjp is calculated from the sum of the momenta of 
the fermions / and /'. If no photon recombination is applied, always the bare momenta 
are taken. In the case of photon recombination the momentum of recombined photons is 
included in the invariant mass as described in Section n~T\ 

For Mh = 170 GeV, where the width is dominated by contributions where both in- 
termediate W bosons are simultaneously resonant, the shape of the curves in Figure ITUl 
can be understood as follows. If one of the fermions resulting from the decay of a res- 
onant W boson emits a photon, the invariant mass Mjj, is reduced, giving rise to an 
enhancement for small invariant masses. Without photon recombination these positive 
corrections are large due to the appearance of logarithms of the small fermion masses. 
As the electron mass is smaller, the corresponding logarithms yield a larger contribution. 
If photon recombination is applied, events are rearranged from small invariant masses 
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Figure 9: Predictions for the partial decay widths for H -^ v^e^ ^"v^j^ and H — *■ e~e'^fi^fi^ 
obtained with the program Hdecay normahzed to the complete lowest-order decay width. 
The corrections shown in Figures [7| and |H1 are included for comparison. 
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to large invariant masses. In this case, the observable is inclusive, i.e. the fermion mass 
logarithms cancel owing to the KLN theorem, and the t'eC"'" and /i~z/^ distributions do 
not differ. The analogous phenomenon has, e.g., been discussed for the related resonance 
processes e+e" -^ WW -^ 41eptons [ El 1101 and 77 ^ WW ^ 4/ [EH]- 

For Mh = 140 GeV, i.e. below the threshold, only one W boson can become on shell. 
Thus, there is still a resonance around M^p ~ Mw, but also an enhancement below an 
invariant mass of Mh ~ ^w ~ 60 GeV, where the other decaying W boson can become 
resonant. Near the resonance at M^j, ~ Mw the corrections look similar to the doubly- 
resonant case discussed for Mh = 170 GeV above. The same redistribution of events from 
higher to lower invariant mass due to FSR happens as explained above. Between My/ and 
Mh — My/ none of the W bosons is resonant, and the contribution to the lowest-order 
width is small. Therefore, owing to the redistribution of events via photon emission the 
relative corrections are large in this regime. Below Mh — M-w, where the other W boson 
can become resonant, qualitatively the same FSR effects are visible as in the vicinity of 
the resonance at My/: apart from a constant positive off-set in the relative corrections, 
events are distributed from the right to the left of the maximum. 

Figure ^2 shows the corresponding invariant-mass distributions for the decay H —>■ 
e~e'^fi~fi^ with Mh = 200 GeV and Mh = 170 GeV. The generic features of the plots are 
similar to the decay into W bosons. For Mh = 200 GeV, i.e. above the ZZ threshold, there 
is a resonance region around Mz, and the corrections become large in the non-collinear- 
safe case. Photon recombination rearranges the events, so that the fermion logarithms 
cancel. For Higgs masses below the ZZ threshold, such as for Mh = 170 GeV, one Z boson 
or the other is resonant for M^j- ~ Mz or Mjj < Mh — Mz, respectively. The shape and 
the large size of the corrections are due to collinear FSR as explained above. In Ref. [ ITHj 
it was pointed out that the kinematical threshold near Mh — Mz where the other Z boson 
can become on shell, which is at Mjf ~ 80 GeV in Figure ^J can be used to verify the 
spin of the Higgs boson. While the rise of the width near this threshold is proportional 
to /? oc ^'[l - (Mz + Mj/)VMg][l - (Mz - Mff-y/M^] for a spin-0 particle, it would be 
proportional to (3^ for a spin-1 particle. Figure ITT] shows that the radiative corrections 
influence the slope at the kinematical threshold. 

Finally, in Figure ^1 we investigate the influence of the contribution (5fsr of higher- 
order FSR to the complete relative correction 6 on the invariant-mass distribution of 
/i~z/^ and fi~fi'^ in the decays H — >■ i>ee~^fi~Ufj_ and H — >■ e~e"'"/i~/i"'". The invariant mass 
is defined via the momenta of the fermions alone, i.e. without photon recombination. If 
photon recombination was applied, the leading logarithmic FSR corrections, as described 
in Section 1131 would vanish completely. Subtracting the 0{a) terms ()4.16|) with ()4.17p 
from (J4.1H) with the structure functions (|4.12j) yields the contribution that is beyond 
0{a). In Figure IT^ the impact of this contribution is studied revealing corrections of up 
to 4% in regions where the lowest-order result is relatively small. Figure ^1 also shows 
the comparison between the structure function with and without the exponentiation of 
the soft-photon parts in ()4.12j) and ()4.14|) . respectively. The difference is beyond 0{a^) 
and turns out to be tiny. 
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Figure 11: Distribution in the invariant mass of the /x /i+ pair (l.h.s.) and relative cor- 
rections to the distributions in the invariant masses of the e~e"'" and fi~fi~^ pairs (r.h.s.) 
in the decay H -^ e-e+/i-/i+ for Mh = 200 GeV and Mh = 170 GeV. 
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Figure 12: Influence of the leading logarithmic terms of FSR on the invariant-mass distri- 
bution oi fi~Un and fi~fi~^ in the decays H — >■ Ve^'^ ^~v^ and H — >■ e~e"'"/i~/i^. The different 
curves correspond to the result with exponentiation, without exponentiation, and to the 
sum of a^ and a^ terms, which are labelled "beyond (9(a)". 

7.5 Angular distributions 

The investigation of angular correlations between the fermionic decay products is an 
essential means of testing the properties of the Higgs boson. In Refs. [ [T21 CSl it was 
demonstrated how the spin of the Higgs boson can be determined by looking at the angle 
between the decay planes of the Z bosons in the decay H — i> ZZ. This angle can be defined 
by 



cos I 



sgn(smi 



(ki2 X ki)(ki2 X kg) 
|ki2 X ki||ki2 X ksl ' 
sgn{ki2 ■ [(ki2 X ki) x (ki2 x kg)]}. 



(7.9) 



where ki2 



= ki + k2. The l.h.s. of Figure IT^ shows the differential decay width for H — > 
as a function of 0' revealing a cos 20' term. As was noticed in Refs. [ ^J E] , 
this term would be proportional to (— cos 20') if the Higgs boson was a pseudo-scalar. 

Note that for non-photonic events the definition of 0' coincides with the definition 
given in Ref. [EHl where (— k34 x kg) with k34 = kg -|- k4 was used instead of (ki2 x kg). 
Explicitly, was defined by 



COS( 



sgn(sm( 



(ki2 X ki)(-kg4 X kg) 
|kl2 X kill -kg4 X kg|' 

sgn{ki2 ■ [(ki2 X ki) x 



-kg4 X kg)]}. 



(7.10) 



However, this definition yields large negative contributions at = 0° and = 180°. As 
was explained in Ref. [EH], this is an effect of the suppressed phase space in the real 
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Figure 13: Distribution in the angle between the Z —>■ I l~^ decay planes in the decay H — > 
e~e~^fi~fi~^ (l.h.s.) and corresponding relative corrections (with photon recombination) 
(r.h.s.) for Mh = 200 GeV. 

corrections. At = 0° and (p = 180° the phase space for photonic events shrinks to the 
configurations where the photon is either soft or lies in the decay plane of the gauge bosons. 
Thus, the negative contributions from the virtual corrections are not fully compensated 
by the real corrections. Using ki2 x ka as in ()7.9|) avoids this suppression and gives rise 
to a smooth dependence of the corrections on as can be seen on the r.h.s. of Figure El 
which shows the relative corrections for and 0' in the decay H — >■ e~e"'"/x~/i"'". Since the 
difference of and 0' is only due to photons, this, again, emphasizes the influence of the 
photon treatment. 

In contrast to the invariant-mass distribution of Figure ^J photon recombination does 
not produce any signiflcant effect for the observables 0, 0'. This is because adding a soft 
or collinear photon to a fermion momentum does not change its direction signiflcantly 
and, thus, has only a small influence on the angles 0, 0'. 

The distribution in the decay angle of the /i~ relative to the corresponding Z boson in 
the decay H — > e~e'^fi~fi~^ is shown in Figure El The angle is deflned in the rest frame of 
the Z boson. Since the Z bosons resulting from Higgs decay are preferably longitudinally 
polarized, the distribution involves a component proportional to sin^ ^z^- . The relative 
corrections which are shown in the plot on the r.h.s. reveal a strong enhancement in 
the forward and backward direction if no recombination is applied. This enhancement is 
due to events where the /i"*" emits a collinear photon and has only a small energy left. 
Since the momentum of the Z boson is deflned via its decay fermions, it has almost the 
same momentum as the /i~. This conflguration is enhanced by collinear logarithms which 
are not compensated by virtual contributions. After applying photon recombination, the 
momentum of the Z boson is deflned via the sum of the fermion and photon momenta. 
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Figure 14: Distribution in the angle between the yU~ and the corresponding Z boson in 
the rest frame of the Z boson (l.h.s.) and corresponding relative corrections with and 
without photon recombination (r.h.s.) in the decay H — > e~e^/i~/i^ for Mh = 170 GeV 
and Mh = 200 GeV. 
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Thus, the fi~ is not necessarily colhnear to the Z boson anymore, and events are rearranged 
to smaller | cos 9z^- \ giving rise to a flatter distribution. 

Next, we consider the distribution in the angle between two fermions. In the case of 
H —>■ WW the angle between the charged fermions can be used to discriminate the Higgs 
signal events from background events [E], because the fermions are emitted preferably in 
the same direction. This can be understood as follows. At leading order, the only non- 
vanishing helicity amplitudes for H —>■ WW are those with equal-helicity W bosons. Since 
W bosons only couple to left-handed particles and due to angular momentum conservation, 
particles (anti-particles) are emitted preferably in the forward direction of transverse 
W bosons with negative (positive) helicity, and anti-particles (particles) in the backward 
direction. As, close to threshold, 2/3 of the W bosons are transverse and as the W bosons 
fly in opposite directions, a particle and an anti-particle of their decay products will be 
emitted preferably in the same direction, resulting in small angles between these particles. 

In the decay H -^ v^^^ ^'^v^ neither the Higgs-boson nor the W-boson momenta can 
be reconstructed from the decay products. The distribution in the angle between the e+ 
and /i~ can, thus, only be studied upon including the Higgs-production process. If the 
Higgs boson was, however, produced without transverse momentum, or if the transverse 
momentum was known, the angle between e^ and ^~ in the plane perpendicular to the 
beam axis could be studied without knowledge of the production process. We define the 
transverse angle between e^ and /i~ in the frame where kH,T = as 

k2 T ■ ka T 

COS0eM,T = Tj TTj f, 

sgn(sin0e/.,T) = sgnje^ ■ (k2,T x kg^x)}, (7.11) 

where kj t are the transverse components of the fermion momenta w.r.t. the unit vector 
e^, which could be identified with the beam direction of a Higgs production process. 

The corresponding distribution, together with the influence of the corrections, is shown 
in Figure IT^ The enhancement for small angles, which was explained above, is transferred 
to the distribution of the transverse angle 0e/i,T- Since the photon recombination does 
not change the direction of the fermions, it does not have any visible effect on the relative 
corrections. 

Finally, we investigate the distribution of the angle between e~ and yU~ in the decay 
H — > e~e^/i~;U"*". We prefer to choose the angle between two fermions with the same 
charge because this constitutes an unambiguous choice in the decay H — > yU~/i+/i~/i+. 
Figure ^1 shows the tendency that the fermions are emitted in opposite directions for 
the same reason as explained above. However, this feature is not as pronounced as in 
H -^ h'ee'^^~h'^, because Z bosons do also couple to right-handed fermions so that one 
Z boson might decay into a left-handed fermion and the other into a right-handed fermion. 
The radiative corrections tend to reduce the enhancement in forward direction and do not 
depend on photon recombination. 

8 Conclusions 

The decays of the Standard Model Higgs boson into four leptons via a W-boson or 
Z-boson pair lead to experimental signatures at the LHC that are both important for the 
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Figure 15: Distribution in the transverse angle between e"*" and /i~ including corrections 
(l.li.s.) and corresponding relative corrections (r.h.s.) with and without applying photon 
recombination in the decay H -^ Uce^fi^u^ for Mh = 140 GeV and Mh = 170 GeV. 
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Figure 16: Distribution in the angle between e~ and /i~ including corrections (l.h.s.) and 
corresponding relative corrections (r.h.s.) with and without applying photon recombina- 
tion in the decay H — > e~e"'"/i~/^"'" for Mr = 200 GeV. 

search for the Higgs boson and for studying its properties. To exploit this possibility a 
Monte Carlo event generator for the decays H -^ WW/ZZ -^ 41eptons is needed that 
properly accounts for the relevant radiative corrections, in order to achieve the necessary 
precision in predictions. Prophecy4f is an event generator dedicated to this task. We 
have shown first results of this generator and described the underlying calculation. 

In detail, we have presented the complete electroweak radiative corrections of 0{a) 
to the decays H -^ 41eptons, supplemented by corrections beyond 0{a) originating from 
heavy-Higgs effects and final-state radiation. The intermediate W- and Z-boson reso- 
nances are treated in the so-called complex-mass scheme, which fully preserves gauge 
invariance and does not employ any type of expansion or on-shell approximation for 
the intermediate gauge-boson resonances. Consequently, the calculation is equally valid 
above, in the vicinity of, and below the WW and ZZ thresholds. 

The corrections to partial decay widths typically amount to some per cent and increase 
with growing Higgs mass Mh, reaching about 8% at Mh ~ 500 GeV. This statement, 
however, applies only if the lowest-order decay widths are already evaluated with the full 
off-shell effects of the intermediate W and Z bosons, in particular near and below the 
WW and ZZ thresholds. The on-shell (narrow-width) approximation for the corrections 
is good within 0.5—1% of the width for Higgs masses sufficiently above the corresponding 
gauge-boson pair threshold, as long as the lowest-order prediction consistently includes 
the off-shell effects of the gauge bosons. For H -^ WW — *> 4/ the narrow-width approxi- 
mation fails by about 10% for Higgs masses that are only 2 GeV above the WW threshold, 
because the instability of the W bosons significantly infiuences the Coulomb singularity 
near threshold. Only a calculation that keeps the full off-shellness of the W and Z bosons 
can describe the threshold regions properly. We have given a simple improved Born ap- 
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proximation for the partial widths that reproduces the full calculation within < 2% for 
Higgs masses below 400 GeV. In this regime our complete calculation has a theoretical 
uncertainty below 1%. For larger Higgs masses we expect that unknown two-loop correc- 
tions that are enhanced by G^M^ deteriorate the accuracy. Finally, for Mh ^ 700 GeV 
it is well known that perturbative predictions become questionable in general. 

For angular distributions, which are important in the verification of the discrete quan- 
tum numbers of the Higgs boson, the corrections are of the order of 5—10% and distort 
the shapes. For invariant-mass distributions of fermion pairs, which are relevant for the 
reconstruction of the gauge bosons, the situation is similar to gauge-boson pair production 
processes such as e"'"e~ — >■ WW — *> 4fermions, i.e. the corrections can reach several tens 
of per cent depending on the treatment of photon radiation. 

In its present version the Monte Carlo event generator Prophecy4f deals with fully 
leptonic final states, a situation most relevant for the LHC The generalization to semi- 
leptonic and hadronic final states, including a proper description of QCD corrections, will 
be described in a forthcoming publication. 

Acknowledgements 

We thank M. Spira for helpful discussions about Hdecay. 

References 

[1] E. W. N. Glover, J. Ohnemus and S. S. D. Willenbrock, Phys. Rev. D 37 (1988) 
3193; 

V. D. Barger, G. Bhattacharya, T. Han and B. A. Kniehl, Phys. Rev. D 43 (1991) 
779; 

V. D. Barger, R. J. N. Phillips and D. Zeppenfeld, Phys. Lett. B 346 (1995) 106 
|hep-ph/9412276|; 

D. L. Rainwater and D. Zeppenfeld, Phys. Rev. D 60 (1999) 113004 [Erratum-ibid. 
D 61 (2000) 099901] hep-ph/9906218|; 

N. Kauer, T. Plehn, D . L. Rainwater and D. Zeppenfeld, Phys. Lett. B 503 (2001) 
113 |hep-ph/001235Tl. 



[2] M. Dittmar and H. K. Dreiner, Phys. Rev. D 55 (1997) 167 |hep-ph/96083T7 



[3] S. Asai et a/., Eur. Phys. J. C 32S2 (2004) 19 |hep-ph/0402254"; 
S. Abdullin et al, Eur. Phys. J. C 39S2 (2005) 41. 

[4] L. Zivkovic, Czech. J. Phys. 54 (2004) A73. 

[5] ATLAS Collaboration, Technical Design Report, CERN-LHCC 99-14 (May 1999); 
CMS Collaboration, Technical Proposal, CERN-LHCC 94-38 (Dec. 1994). 

[6] K. A. Assamagan et al. [Higgs Working Group Collaboration], proceedings of 
the 3rd Les Houches Workshop: "Physics at TeV Colliders", Les Houches, 2003, 



hep-ph/040615^ 



37 



[7] C. Buttar et al. [SMH Working Group Collaboration] QCD, EW, and Higgs working 
group report of the workshop "Physics at TeV Colliders" , Les Houches, May 2005, 
|hep-ph/0604120| 

[8] J. A. Aguilar-Saavedra et al., TESLA Technical Design Report Part III: Physics at 
an e"'"e~ Linear Collider, hep-ph/0106 315| 



[9] T. Abe et al. [American Linear Collider Working Group Collaboration], in Proc. of 
the APS/DPF/DPB Summer Study on the Future of Particle Physics (Snowmass 
2001) ed. R. Davidson and C. Quigg, SLAC-R-570, Resource hook for Snowmass 
2001, Ihep-ex/O lOeOSS, |hep-ex /0106056 , hep-ex/0106057[ |hep-ex /0106058 . 

[10] K. Abe et al. [AGFA Linear Collider Working Group Collaboration], AGFA Linear 
Collider Working Group report, |hep-ph/0109166|. 

[11] N. Meyer and K. Desch, Eur. Phys. J. C 35 (2004) 171. 

[12] C. A. Nelson, Phys. Rev. D 37 (1988) 1220; 

A. Soni and R. M. Xu, Phys. Rev. D 48 (1993) 5259 Lhep-ph/9301225j ; 

D. Chang, W. Y. Keung and I. Phillips, Phys. Rev. D 48 (1993) 3225 

[hep-ph/9303226|; 

A. Skjold and P. Osland, Phys. Lett. B 311 (1993) 261 |hep-ph/9303294"; 

V.D. Barger, K.M. Cheung, A. Djouadi, B.A. Kniehl and P.M. Zerwas, Phys. Rev. 

D 49 (1994) 79 hep-ph/9306270j ; 

T. Arens and L.'m. Sehgal, Z. Phys. C 66 (1995) 89 |hep-p h/9409396 ; 

C. P . Buszello, I . Fleck, P. Marquard and J. J. van der Bij, Eur. Phys. J. C 32 (2004) 

209 |hep-ph/0212396| . 

[13] S. Y. Choi, D. J. Miller, M. M. Miihlleitner and P. M. Zerwas, Phys. Lett. B 553 



(2003) 61 hep-ph/0210077|. 



[14] G. Pocsik and T. Torma, Z. Phys. C 6 (1980) 1; 
T. G. Rizzo, Phys. Rev. D 22 (1980) 722; 
W. Y. Keung and W. J. Marciano, Phys. Rev. D 30 (1984) 248. 

[15] R. N. Cahn, Rept. Prog. Phys. 52 (1989) 389; 
B. A. Kniehl, Phys. Lett. B 244 (1990) 537; 

A. Grau, G. Panchieri and R. J. N. Phillips, Phys. Lett. B 251 (1990) 293; 
E. Gross, G. Wolf and B. A. Kniehl, Z. Phys. C 63 (1994) 417 [Erratum-ibid. C 66 
(1995) 321] ihep-ph/9404220j . 

[16] A. Djouadi, |h^-ph/0503172| 

[17] A. Djouadi, J. Kalinowski and M. Spira, Comput. Phys. Commun. 108 (1998) 56 
[hep-ph/97 04448i . 

[18] J. Fleischer and F. Jegerlehner, Phys. Rev. D 23 (1981) 2001; 
B.A. Kniehl, Nucl. Phys. B 352 (1991) 1; 

D.Y. Bardin, P.K. Khristova and B.M. Vilensky, Sov. J. Nucl. Phys. 54 (1991) 833 
[Yad. Fiz. 54 (1991) 1366]. 

38 



[19] B.A. Kniehl, Nucl. Phys. B 357 (1991) 439. 

[20] B. A. Kniehl and M. Spira, Z. Phys. C 69 (1995) 77 [hep-ph/9505225l;^ 

B. A. Kniehl and M. Steinhauser, Phys. Lett. B 365 (1996) 297 hep-ph/95073821; 
A. Djouadi, P. Gambino and B. A. Kniehl, NucL Phys. B 523 (1998) 17 
[hep-ph/9712330|. 



[21] B. A. Kniehl and M. Steinhauser, Nucl. Phys. B 454 (1995) 485 |hep-ph/9508"24T 
[22] A. Ghinculov, Nucl. Phys. B 455 (1995) 21 |hep-ph/ 9507240|. 



[23] A. Frink, B. A. Kniehl, D. Kreimer and K. Riesselmann, Phys. Rev. D 54 (1996) 
4548 |hep-ph/ 9606310|. 

[24] C M. Carloni Calame, M. Moretti, G. Montagna, O. Nicrosini, F. Piccinini and 
A. D. Polosa, Nucl. Phys. Proc. Suppl. 157 (2006) 73 hep-ph/0604033 1 . 

[25] A. Bredenstein, A. Denner, S. Dittmaier and M.M. Weber, Nucl. Phys. Proc. Suppl. 
157 (2006) 63 | he p-ph/0604033, . 

[26] G. Belanger, F. Boudjema, J. Fujimoto, T. Ishikawa, T. Kaneko, K. Kato and 
Y. Shimizu, Nucl. Phys. Proc. Suppl. 116 (2003) 353 f hep-ph/0211268| and Phys. 
Lett. B 559 (2003) 252 |h ep-ph/0212261, . 

[27] A. Denner, S. Dittmaier, M. Roth and M. M. Weber, Phys. Lett. B 560 (2003) 196 
^hep-ph/030118 9| and Nucl. Phys. B 660 (2003) 289 ^hep-ph/0302198, . 

[28] A. Denner, S. Dittmaier, M. Roth and L. H. Wieders, Phys. Lett. B 612, 223 (2005) 



|hep-ph/0502063i and Nucl. Phys. B 724 (2005) 247 |hep-ph/0505042| . 

[29] A. Denner, S. Dit tmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 560 (1999) 33 
|hep-ph/9904472 . 



[30] A. Denner and S. Dittmaier, Nucl. Phys. B 658 (2003) 175 jhep-ph/02 12259 



[31] A. Denner and S. Dittmaier, Nucl. Phys. B 734 (2006) 62 ^hep-ph/0509141 



[32] S. Dittmaier, Nucl. Phys. B 565 (2000) 69 hep-ph/9904440 



[33] A. Bre denstein, S. Dittmaier and M. Roth, Eur. Phys. J. C 44 (2005) 27 
| hep-ph /0506005|. 

[34] A. Denner, Fortsch. Phys. 41 (1993) 307. 

[35] A. Denner, S. Dittmaier, and G. Weiglein Nucl. Phys. B 440 (1995) 95 
^hep-ph/9410338| . 



[36] S. Dittmaier, Phys. Rev. D 59 (1999) 016007 |hep-ph/9805445 
[37] W. F. L. HoUik, Fortsch. Phys. 38 (1990) 165. 



39 



[38] J. Kiiblbeck, M. Bohm and A. Denner, Comput. Pliys. Coinmun. 60 (1990) 165; 
H. Eck and J. Kiiblbeck, Guide to FeynArts 1.0, University of Wiirzburg, 1992. 



T. Hahn, Comput. Phys. Commun. 140 (2001) 418 hep-ph/0012260 



[40] T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153 
|hep-ph/9807565|; 
T. Hahn, Nucl. Phys. Proc. Suppl. 89 (2000) 231 |hep-ph/0005029| . 



[41] M. W. Griinewald et ai, in Reports of the Working Groups on Precision Galculations 
for LEP2 Physics, eds. S. Jadach, G. Passarino and R. Pittau (CERN 2000-009, 
Geneva, 2000), p. 1 |hep-ph/0005309^ . 

[42] G. 't Hooft and M. Veltman, Nucl. Phys. B 153 (1979) 365. 

[43] W. Beenakker and A. Denner, Nucl. Phys. B 338 (1990) 349. 

[44] A. Denner, U. Nierste and R. Scharf, Nucl. Phys. B 367 (1991) 637. 

[45] G. Passarino and M. Veltman, Nucl. Phys. B 160 (1979) 151. 

[46] A. Sirlin, Phys. Rev. D 22 (1980) 971. 



[47] T. Stelzer and W.F. Long, Comput. Phys. Commun. 81 (1994) 357 | hep-ph/9401258| . 

[48] T. Kinoshita, J. Math. Phys. 3 (1962) 650; 

T. D. Lee and M. Nauenberg, Phys. Rev. 133 (1964) B1549. 

[49] E. A. Kuraev and V. S. Fadin, Sov. J. Nucl. Phys. 41 (1985) 466 [Yad. Fiz. 41 (1985) 
733]; 

G. Altarelli and G. Martinelli, in Physics at LEP, eds. J. Ellis and R. Peccei, (CERN 
86-02, Geneva, 1986), Vol. 1, p. 47; 

O. Nicrosini and L. Trentadue, Phys. Lett. B 196 (1987) 551 and Z. Phys. C 39 
(1988) 479; 

F. A. Berends, W. L. van Neerven and G. J. H. Burgers, Nucl. Phys. B 297 (1988) 
429 [Erratum-ibid. B 304 (1988) 921]. 

[50] W. Beenakker et al, in Physics at LEP2 (Report CE RN 96-01, Genev a, 1996), eds. 
A. Altarelli, T. Sjostrand and F. Zwirner, Vol. 1, p. 3, hep-ph/9602351[ 



[51] F. A. Berends, R. Pittau and R. Kleiss, Nucl. Phys. B 424 (1994) 308 
|hep-ph/9404313| and Comput. Phys. Commun. 85 (1995) 437 hep-ph/9409326]; 
F. A. Berends, P. H. Daverveldt and R. Kleiss, Nucl. Phys. B 253 (1985) 441; 
J. Hilgart, R. Kleiss and F. Le Diberder, Comput. Phys. Commun. 75 (1993) 191. 

[52] A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Comput. Phys. Commun. 153, 
462 (2003) |hep-ph/0209330 . 



[53] A. Bredenstein, S. Dittmaier and M. Roth, Eur. Phys. J. C 36 (2004) 341 



hep-ph/0405169| . 

40 



[54] G.P. Lepage, J. Comput. Pliys. 27 (1978) 192. 

[55] V. S. Fadin, V. A. Khoze and A. D. Martin, Phys. Lett. B 311 (1993) 311; 
D. Y. Bardin, W. Beenakker and A. Denner, Phys. Lett. B 317 (1993) 213. 

[56] S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B 592 (2004) 1. 



[57] F. Jegerlehner, |hep-ph/0105283, 



[58] D. Y. Bardin, A. Leike, T. Riemann and M. Sachwitz, Phys. Lett. B 206 (1988) 539; 
D. Wackeroth and W. Hollik, Phys. Rev. D 55 (1997) 6788 hep-ph/9606398j ; 
W. Beenakker et al, Nucl. Phys. B 500 (1997) 255 |hep-pli/9612260r 



A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 587 (2000) 67 
|he p-ph/0006307j . 

[60] W. Beenakker, F. A. Berends and A. P. Chapovsky, Nucl. Phys. B 548 (1999) 3 
^hep-ph/9811481j . 



41 



